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Searches for extrasolar planets have uncovered an astonishing diversity of 
planetary systems, yet the frequency of solar system analogs remains unknown. 
The gravitational microlensing planet search method is potentially sensitive to 
multiple-planet systems containing analogs of all the solar system planets ex- 
cept Mercury. We report the detection of a multiple-planet system with mi- 
crolensing. We identify two planets with masses of ~ 0.71 and ~ 0.27 times 
the mass of Jupiter and orbital separations of ~ 2.3 and ~ 4.6 astronomi- 
cal units orbiting a primary star of mass ~ 0.50 solar masses at a distance of 
~ 1.5 kiloparsecs. This system resembles a scaled version of our solar system 
in that the mass ratio, separation ratio, and equilibrium temperatures of the 
planets are similar to those of Jupiter and Saturn. These planets could not 
have been detected with other techniques; their discovery from only six con- 
firmed microlensing planet detections suggests that solar system analogs may 
be common. 

Nearly 250 extrasolar planets (i ) have been discovered by measuring a variety of effects: 
reflex motion of the host star using pulsar timing or precision Doppler measurements d2]|21El); 
periodic dimming of the parent star as the planet transits in front ([5]|6|); and planet-induced 
perturbations to microlensing light curves in which the host star acts as the primary gravitational 
lens ([71IS1I91 I701I77]) . These detections have uncovered an enormous range of planetary properties, 
indicating that planetary systems very unlike our own are common throughout the Galaxy (fi2]) . 

To date, ~ 25 multiple-planet systems have been detected (1731) . all but one (2) using the 
Doppler method. Because Doppler surveys must monitor the host star's reflex motion over the 
planet's orbital period, they are limited by the finite duration as well as the sensitivity of the 
measurements. Hence, they are only just now becoming sensitive to Jupiter analogs and are 
not yet sensitive to Saturn analogs (nor, ipso facto, Jupiter/Saturn systems). Thus, all multiple- 
planet systems discovered so far are very dissimilar from our own, and the frequency of solar 
system analogs remains unknown. 

Because microlensing relies on the direct perturbation of light from distant stars by the 
gravitational field of the planet, it is 'instantaneously' able to detect planets without requiring 
observations over a full orbit. For a primary star of mass M, microlensing sensitivity peaks 
for planets in the range ~ [1 — 5](M/0.3M Q ) 1 / 2 astronomical units (AU) For solar-mass 
stars, this is exactly the range of the solar system gas giants so microlensing provides a method 
to probe solar system analogs ([7411751) . 
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As pointed out by Griest & Safizadeh (16), the very rare class of high-magnification (> 100) 
microlensing events provides an extremely sensitive method of detecting planets. Near the peak 
of high-magnification events, the two images created by the primary star are highly magnified 
and distorted, and form a complete or nearly complete Einstein ring. A planetary companion 
to the primary star lying reasonably near the Einstein ring will distort the symmetry of the 
ring. As the host passes very close to the source line-of- sight, the images sweep around the 
Einstein ring, thus probing this distortion. Although the total number of high-magnification 
events is small, the instantaneous chance of detection in each is much higher than for the more 
common low-magnification events. Equally important, the interval of high-sensitivity (i.e., 
high-magnification) is predictable from the evolution of the light curve (|76]|77l[7^]|79l) . This 
permits concentration of scarce observing resources on these events. Furthermore, the high- 
magnification makes it possible to acquire high signal-to-noise ratio photometry of the peak of 
the events using relatively small-aperture (and so plentiful) telescopes. As a result, four ([911771) 
of the six planets (IffllTOl) discovered to date in microlensing events were in high-magnification 
events. 

Almost immediately after Griest and Safizadeh (16) pointed out the sensitivity of high- 
magnification events, Gaudi et al. ([201) derived an important corollary. Because planets in the 
neighborhood of the Einstein ring are revealed with near unit probability in high-magnification 
events, multiple-planet systems lying in this region will be revealed with almost the same prob- 
ability. 

The Optical Gravitational Lens Experiment (OGLE) ([27]) and Microlensing Observations 
in Astrophysics (MOA) ([791) collaborations together alert ~ 700 ongoing microlensing events 
per year. Two collaborations, a joint venture of the Probing Lensing Anomalies NETwork 
(PLANET) ([221) and RoboNet (|25l) collaborations, and the Microlensing Follow-Up Network 
(//FUN) (24), then monitor a subset of these alerts to search for planets. //FUN focuses almost 
entirely on high-magnification events, including two events originally alerted by OGLE that 
proved to have a Jupiter-mass (9) and a Neptune-mass (11) planet, respectively. Here, we report 
on the detection of a multi-planet system using this approach. 

On 28 March 2006 (HJD~ 3822), the OGLE Early Warning System (EWS) <|27J) announced 
OGLE-2006-BLG-109 as a non-standard microlensing event possibly indicative of a planet. 
This immediately triggered followup observations by //FUN and RoboNet, which gained inten- 
sity as the event approached high-magnification. On 5 April, the event underwent a deviation 
from the single-lens form indicative of a binary lens. Within 12 hours of this deviation, a pre- 
liminary model indicated a jovian-class planet, which was predicted to generate an additional 
peak on 8 April. The 8 April peak occurred as predicted, but in the meantime, there was an 
additional peak on 5/6 April, which turned out to be due to a second Jovian-class planet. 

Figure 1 shows the data from 11 observatories, including 7 from //FUN [the Auckland 
0.35m and Farm Cove 0.25m in New Zealand (clear filter), the Wise lm in Israel (clear), the 
CTIO/SMARTS 1.3m in Chile (/-band and if -band), the Areo8 0.3m in New Mexico operated 
by the Campo Catino Astronomical Observatory (clear), and the MDM 2.4m (f-band) and Mt. 
Lemmon 1.0m (f-band) in Arizona], the OGLE Warsaw 1.3m (f-band) in Chile, the MOA Mt. 
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John 0.6m (/-band) in New Zealand, the PLANET/Canopus lm (/-band) in Tasmania, and the 
RoboNet/Liverpool 2m (/?-band) in the Canary Islands. There are a total of 2642 data points. 
In addition, there are 29 V-band data points from OGLE and CTIO/SMARTS that we use to 
determine the source color. 

The qualitative character of the event can be read directly from the light curve, primar- 
ily from the five distinctive features shown in Figure 2. Consider the first three features: the 
low-amplitude anomaly (OGLE, HJD~ 3823) that triggered the OGLE EWS alert, the gentle 
"shoulder" during the first rise (MOA, HJD~ 3830), and the first peak (Auckland, HJD~ 3831). 
Together, these can only be produced by, respectively, passage close to or over a weak cusp, en- 
trance into a weak caustic, and exit from a strong caustic. (The magnification diverges when 
a point source crosses a closed concave caustic curve, where additional images are created on 
entry or destroyed on exit of the enclosed area. Caustics are strong or weak depending on 
the brightness of these images. The concave curves meet at cusps that produce sharp spikes 
of magnification when crossed by the source.) Such a sequence requires a topology similar 
to the one shown in the inset to Figure 1. The specific strengths of each feature require the 
specific caustic topology shown. In particular, the narrow mouth of the caustic toward the bot- 
tom generates a very strong caustic. This was essentially the argument used to predict the fifth 
feature (OGLE/MDM/Lemmon/Auckland/FarmCove/Tasmania, HJD~ 3834), corresponding 
to a moderately strong cusp passage (Fig. 1). The size and strength of this caustic imply a 
jovian-class planet lying very close to the Einstein ring, although detailed modeling is required 
to derive the precise planet/star mass ratio. The fourth feature (Wise/OGLE, HJD~ 3831.5) 
cannot be explained by considering the caustic generated by this jovian-class planet alone. This 
feature occurs near the time when the source approaches closest to the center-of-mass of the 
planet/star system; this is exactly the time at which the central-caustic bumps due to additional 
planets are expected to occur (|20l) . The inset in Figure 1 highlights the additional caustic fea- 
ture due to a second planet that is required to explain this bump. This caustic feature is smaller 
than the main caustic, which implies that the planet, also of jovian class, lies farther from the 
Einstein ring, so it is subject to the standard (|7o1) b *-* b^ 1 degeneracy, where b is the planet-star 
projected separation in units of the Einstein radius. A detailed analysis shows the mass is three 
times as great as that of the first planet and that the b < 1 solution is favored by A% 2 = 11.4. 
We label these planets OGLE-2006-BLG-109Lc and OGLE-2006-BLG-109Lb, respectively. 
Although the caustics of the individual planets do interact to form a single caustic curve, their 
effects are nevertheless mostly independent (1251 I7ffll2"6l) . so the parts of the caustic associated 
with the individual planets can be identified, as shown in Figure 1 . Modeling the light curve in 
detail with a three-body lens yields, m^jM = 1.35 x 10~ 3 , m c /m b = 0.36 for the mass ratio of 
the planets and their host, very similar to rrij/M & = 0.96 x 10" 3 and m s /rrij = 0.30 for Jupiter, 
Saturn, and the Sun. The ratio of projected separations r± t b/r± tC = 0.60 is also very similar to 
the Jupiter/Saturn value of aj/a s = 0.55. 

Two subgroups of authors conducted independent searches for alternative solutions. Both 
found that no single-planet solution is consistent with the light-curve topology. We also succes- 
sively eliminated each of the five features to see whether the remaining four features could be fit 
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by a single planet. We found that only the elimination of the fourth feature produced successful 
single-planet models. By contrast, similar procedures in other events (|77ll27l) led to many in- 
dependent solutions. OGLE-2006-BLG-109 differs from these in that it has five well-covered 
features. 

Several higher-order effects are apparent in this event that permit us to extract much more 
detailed information about the system from the light curve. We only briefly sketch these here. 
For over 95% of microlensing events observed from the ground, the lens parameters are deter- 
mined only relative to the angular Einstein radius 9 E , whose absolute scale remains unknown. 
Here 6 E = ^AGM/c 2 D, where M is the mass of the lens, l/D = l/D l - 1 /D s , and D x and D s 
are the distances to the lens and the source, respectively. However, in this event the effect of the 
finite size of the source star during caustic exit allows us to measure the source radius relative 
to the Einstein radius, p = O^/Oe (|2^1) . From the source color and flux we can determine its 
angular size 9*, and thus 9 E (|2?)) . 

The acceleration of Earth in its orbit about the Sun induces subtle distortions on the light 
curve called microlens parallax, which yields the physical size of the Einstein radius projected 
onto the observer plane, f E = 9 E D (|29l) . This is usually measured only in the roughly 3% of 
events that are extremely long, but this event happens to be long and so displays clear distortions 
arising from parallax. 

Combining these two measures of the Einstein radius allows us to triangulate the event and 
so determine the host star distance, D t = l/(9 E /f E + 1/D s ), and mass, M = (c 2 /AG)f E 9 E . We 
assume D s = 8 kpc, although our results are insensitive to this assumption. From a preliminary 
analysis we infer D\ ~ 1.5 kpc and M ~ 0.5 M . Based on high-resolution Keck AO H- 
band images, we detect light from the lens and infer its magnitude to be H = 17.17 ± 0.25, 
consistent with the mass estimate from the light curve. We subsequently incorporate the lens 
flux constraint in the light curve analysis, which allows us to derive more precise estimates of 
Di = 1.49 ± 0.13 kpc and M = 0.50 ± 0.05 M . The planet masses are m b = 0.71 ± 0.08 and 
m c = 0.27 ± 0.03 times the mass of Jupiter. 

Finally, we also detect the orbital motion of the outer planet; this motion both rotates and 
changes the shape of the larger caustic shown in the top inset to Figure 1. We are able to con- 
strain the two components of the projected velocity of the planet relative to the primary star. 
Together with the estimate of the stellar mass, they completely determine the outer planet's or- 
bit (including inclination) under the assumption that it is circular, up to a two-fold degeneracy. 
The solution presented here is marginally favored by the data at A% 2 = 4.8 via the effect of 
the planet's acceleration on the light curve. Thus we can estimate the full (three-dimensional) 
separation of planet c (again assuming a circular orbit), and also of planet b (assuming a copla- 
nar and circular orbit). We find a& = 2.3 ± 0.2 AU and a c = 4.6 ± 0.5 AU. A more refined 
estimate of these parameters and their uncertainties will require a detailed analysis including 
the combined effects of finite sources, parallax, and orbital motion of the planets. The results 
of this analysis will be presented elsewhere (Bennett et al., in preparation). 

The OGLE-2006-BLG-109L planetary system bears a remarkable similarity to our own 
solar system. Although the primary mass is only half solar, the mass ratio of the two planets 
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(0.37) and separation ratio (0.50) are similar to those of Jupiter and Saturn. We infer their 
equilibrium temperatures to be T eq ~ 82 ± 12 K and T eq ~ 59 ± 7 K, ~ 30% smaller than those 
of Jupiter and Saturn. 

Before the detection of extrasolar planets, planet formation theories generally predicted that 
other systems should resemble our solar system. In the core- accretion paradigm, the most mas- 
sive giant planet forms at the 'snow line,' the point in the protoplanetary disk exterior to which 
ices are stable. Immediately beyond the snow line, the surface density of solids is highest and 
the dynamical time is the shortest, and therefore the timescale for planet formation is the short- 
est. Beyond the snow line, the formation timescale increases with distance from the host star. 
Thus in this 'classical' picture of planet formation, one would expect planet mass to decrease 
with increasing distance beyond the snow line, as is observed in our solar system (30). The dis- 
covery of a population of massive planets well interior to the snow line demonstrated that this 
picture of planet formation is incomplete, and considerable inward migration of planets must 
occur (TJ7T) . Nevertheless, this classical picture may still be applicable to our solar system and 
some fraction of other systems as well. The OGLE-2006-BLG-109L planetary system repre- 
sents just such a 'scaled version' of our own solar system, with a less-massive host. This system 
preserves the mass-distance correlation in our solar system, and the scaling with primary mass 
is consistent with the core- accretion paradigm in which giant planets that form around lower- 
mass stars are expected to be less massive but form in regions of the protoplanetary disk with 
similar equilibrium temperatures and are therefore closer to their parent star (|32)) . 

The majority of the ~ 25 known multi-planet systems are quite dissimilar to the OGLE- 
2006-BLG-109L system and to our own solar system. Many of these systems have the very 
close-in massive planets indicative of large-scale planetary migration, or they have a 'normal 
hierarchy', in which the masses of the giant planets increase with distance from the parent star. 
There are two multi-planet systems with properties roughly similar to those of OGLE-2006- 
BLG-109L. The 47 UMa and 14 Her systems each contain a giant planet at a semimajor axis 
of ~ 3 AU and a second, less massive giant planet at a separation of ~ 7 AU (TJ5l) . However, 
because of their higher-mass primaries, the equilibrium temperatures of these planets are con- 
siderably higher than those of OGLE-2006-BLG-109L or Jupiter and Saturn, so these systems 
cannot be considered close analogs of our solar system. 

OGLE-2006-BLG-109Lb and OGLE-2006-BLG-109Lc are the fifth and sixth planets to be 
detected by microlensing. Although, given the detection of planet c, the a priori probability 
of detecting planet b in this event was high, it was not unity. Furthermore, only two other 
jovian-mass planets have been detected by microlensing (|#]|9]), and neither event had substan- 
tial sensitivity to multiple planets. These facts may indicate that the stars being probed by 
microlensing that host jovian-mass companions are also likely to host additional giant planets. 
If the OGLE-2006-BLG-109L planetary system is typical, these systems may have properties 
similar to our solar system. Regardless, the detection of the OGLE-2006-BLG-109L planetary 
system demonstrates that microlensing surveys will be able to constrain the frequency of solar 
system analogs throughout the Galaxy. 
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Figure 1: Data and best-fit model of the OGLE-2006-BLG-109Lb,c two-planet system. The 
data have been binned for clarity, although the fitting procedures used the unbinned data. Data 
from each different observatory /filter combination (as indicated by the color scheme) have been 
aligned using a linear fit to the magnification, which introduces negligible uncertainties. Only 
data near the peak of the event are shown (the unlensed magnitude is / = 16.42). Panel A: The 
source trajectory through the caustic created by the two-planet system is shown as the dark grey 
curve with the arrow indicating the direction of motion. The horizontal line shows an angular 
scale of 0.01 6e, or ~ 15 fias. The shape and orientation of the caustic due to both planets at the 
peak of the event is shown by the black curve. The five light-curve features detailed in Fig. [2] 
are caused by the source crossing or approaching the caustic; the approximate locations of the 
features are labeled with numbers. The majority of the caustic (in black) is due to only the outer 
(Saturn-analog) planet; this portion of the caustic explains four of the five features. The portion 
arising from the second (Jupiter-analog) planet is highlighted in red. This additional cusp in the 
caustic is required to explain the fourth feature in the light curve; as such, the fourth feature 
signals the presence of a second (Jupiter-analog) planet. Because of the orbital motion of the 
Saturn-analog planet, the shape and orientation of the caustic changes over the course of the 
event. The light grey curves show the caustic at the time of features 1 and 5. Panel B: A zoom 
of the source trajectory and caustic near the times of the second, third, and fourth features. The 
circle shows the size of the source. 10 
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Figure 2: Five features of light curve from Fig. [Q which determine planetary geometry. A) 
Feature 1 : weak cusp crossing; B) Feature 2: weak caustic entrance; C) Feature 3: strong caustic 
exit; D) Feature 4: strong cusp approach; E) Feature 5: moderate cusp approach. Features 1, 2, 
3, and 5, are explained by the black portion of the caustic seen in in Fig.[l]A Feature 4 requires 
an additional cusp in the caustic, which is shown as the red curve. Data have been binned for 
clarity. 
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